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Abstract. We propose a very long baseline atom interferometer test of Einstein’s 
equivalence principle (EEP) with ytterbium and rubidium extending over 10 m of free 
fall. In view of existing parametrizations of EEP violations, this choice of test masses 
significantly broadens the scope of atom interferometric EEP tests with respect to other 
performed or proposed tests by comparing two elements with high atomic numbers. In 
a first step, our experimental scheme will allow reaching an accuracy in the Eotvos ratio 
of 7 • 10 -13 . This achievement will constrain violation scenarios beyond our present 
knowledge and will represent an important milestone for exploring a variety of schemes 
for further improvements of the tests as outlined in the paper. We will discuss the 
technical realisation in the new infrastructure of the Hanover Institute of Technology 
(HITec) and give a short overview of the requirements to reach this accuracy. The 
experiment will demonstrate a variety of techniques which will be employed in future 
tests of EEP, high accuracy gravimetry and gravity-gradiometry. It includes operation 
of a force sensitive atom interferometer with an alkaline earth like element in free fall, 
beam splitting over macroscopic distances and novel source concepts. 
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Einstein’s equivalence principle (EEP) is at the core of onr understanding of gravitation 
and is among the most important postulates of modern physics. It is under constant 
scrutiny since a violation of any of its pillars would lead to new physics beyond general 
relativity (GR) and would mark an important milestone in the search for a theory of 
everything (TOE). The EEP is comprised of three separate postulates: the Universality 
of Free Fall (UFF), Local Lorentz Invariance (LLI) and Local Position Invariance (LPI). 
Free fall experiments, as the one described in this letter, test the UFF by comparing the 
accelerations of two bodies of different internal structure and mass in a gravitational 
held. This inertial and gravitational mass equality is also known as the weak equivalence 
principle (WEP). To quantify a possible violation of the UFF it is common to normalise 
the acceleration difference between two test masses to the average local gravitational 
acceleration. This parametrization leads to the Eotvos ratio defined by 


Va,b 


2 9a ~ 9b 
9a + 9 b' 


with gA,B being the gravitational acceleration of test masses A and B respectively. The 
most straightforward way to do such a test is to directly measure the acceleration of two 
bodies in the same gravitational held. This class of tests is called Galilean and the most 
accurate to date was performed by comparing uranium and copper at a level of 1CT 10 [1J. 
The most accurate tests of the UFF were performed by the lunar laser ranging project 
(LLR), measuring the free fall of the moon and the earth in the gravitational held of the 
solar system. Since the UFF is a statement about the acting forces, not only Galilean 
type free fall experiments are performed to test it, but also force balance experiments 
with torsion balances. Torsion balances and LLR constrain possible violations of UFF to 
less than 1CU 13 in Eotvos ratio M- No violation was found so far. Future experiments 
with classical bodies are striving towards spaceborne platforms, to reduce the influence 
of external error source and allow measurements far beyond current state of the art m- 
The use of atom interferometry broadens the held of test masses and allows an operation 
in the quantum regime. As such it is a complementary method to experiments with 
macroscopic bodies and will test aspects formerly inaccessible, such as violations linked 
to the coherence length of the test mass jHj, the possibility to employ cold atoms 
as accelerometers and clocks, and the possibility of spin-polarisation [7j. A hrst 
measurement was performed by a device measuring gravity with a fountain of cold 
caesium atoms and comparing their fall rates to a commercial falling corner cube 
gravimeter at a level of 7T0~ 9 [8j. More recent experiments demonstrate tests of the UFF 
by using atom interferometry with two different quantum objects within the same device 
but do not yet reach the same precision. They are in part relying on two isotopes of the 
same species mm® but also on isotopes of two different elements m ■ Especially tests 
with two isotopes want to take benefit from similarities for large noise suppression factors 
intrinsically arising from the measurements arrangement. New experiments of both 
types are proposed to exceed the limits of current sensitivities, either on ground [121 US] 
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or in micro-gravity environments mm, including the STE-QUEST space mission ng. 
To employ this variety of test candidates in a precision experiment, a crucial point is the 
ability to trap both of the species not only simultaneously but rather in the same trap 
to have a well defined overlap of their initial positions and velocities. In this respect we 
propose quantum degenerate mixtures of rubidium and ytterbium for testing the UFF 
in a large scale device on ground. 

In this paper we discuss the unique features of these mixtures that make them an ideal 
choice as test masses by calculating their violation parameters and comparing them to 
the ones used in other experiments and recent proposals. Focusing on the miscibility of 
different isotopes of these two elements, we will give a description on the source setup 
we aim for. Besides this description, we present possible scenarios for performing a UFF 
test with Bragg-type beam splitters. Along this we analyze noise contributions to the 
measured signal and estimate the performance of a test of the UFF to be 7 • 10 -13 in 
the Eotvos ratio. 

2. Choice of test pairs 

As already mentioned the common way of quantifying an experiment testing the UFF 
is the Eotvos ratio, which scales a measured differential acceleration to the strength of 
the local gravitational field comparing any abnormal composition based forces to the 
composition independent force. While this is a reasonable way to quantify the result of 
the performed measurement it does not take into account the specific kind of composition 
dependence in question. By just using the Eotvos parameter as a tool for comparing two 
tests, an experiment with two spin polarized samples of the same isotope would not be 
treated different than a comparison between hydrogen and anti-hydrogen as proposed 
in HZI while being fundamentally different. Taking the specific composition difference 
into account is part of the interpretation of the data and is strongly dependent on the 
model used to assess a possible violation theory. The use of extended wave functions 
for testing UFF opens the path to formerly unexplored theoretical models which are 
probing the quantum nature of matter and its interaction with space time [S]. While 
this is a vast field of study, we will focus on models which allow us a comparison 
to classical experiments. Specifically we asses the dilaton scenario [18] and a scenario- 
independent scaling approach based on the standard model extension (SME) [115]. Atom 
interferometry can provide several new aspects different with respect to classical test 
masses as the test masses are of high isotopic purity and the choices of test masses can 
be extended beyond non-magnetic, conducting solids which are typically used in torsion 
balances. 

According to the dilation model [18] a violation may be caused by forces acting 
differently on neutron and proton number. With the introduced effective charges Qab 
and Q'2 b calculated from the composition of a test particle a measurement of the Eotvos 
ratio set bounds to the parameters D\ and D 2 according to the formula 

??a,b = ^i(^Qa,b) + D 2 (AQ^ B ). 


( 1 ) 
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A similar kind of parameterization can be given for the standard model extension [19] 

hA,B — A/_ n + A f +n + A f_ n + A f +n (2) 

with the defined violation parameters for matter and anti-matter linked to neutron 
excess and total baryon number 

A f-n = f 0l+P -^ e+P ~ n - f^ B+P -nf3 e+P ~ n 

A f+ n = f ne+p+n j3 e+P+n ~ f ne+p+n (3 e+P+n 

/ o \ 

A f — f _, _ - /Qe+p-n _ f _ oe+p—n ^ ' 

L\J_ n ~ J p e +P~ n P J 0^+P~ n P 

A f +n = fp l+P+ nF +P+n - f^ffi****' 

In both models larger absolute differences in the sensitivity factors of the employed 
test mass pair give rise to a larger signal in case of a violation of the UFF. Vice 
versa, an experimental determination of the Eotvos ratio for such a test mass choice 
better constrains the existence of violations than tests performed with lower sensitivity 
factors for the same accuracy. Moreover, different test mass pairs probe different linear 
combinations of suspected violations linked to the neutron excess and the total baryon 
number of the test masses. In order to unambiguously determine the origin of a violation, 
a minimum of two test mass pairs needs to be employed. Interestingly, as shown in 
Ref. [20], even a test performed at a lower accuracy as compared to state of the art tests 
can further constrain possible violations, when the used test masses are significantly 
different to previously utilized ones. The sensitivity factors for different choices of test 
pairs are presented in table [l] For example, in comparison to Be-Ti the combination 
of ytterbium and rubidium isotopes is a factor of 2 more sensitive to baryon number 
related violations and even three orders of magnitude more sensitive in the parameter 

A f- n . 

3. Atom interferometry in a 10 m atomic fountain 

The inertial sensitive interferometry with cold rubidium clouds is well covered by state- 
of-the-art experiments for measuring gravity EU 127!, gravity gradients [2BJ and rota¬ 
tions ra as well as for measuring fundamental constants |30j . Similarly laser-cooled 
ytterbium is by now very successfully utilized in optical clocks, especially optical lat¬ 
tice clocks EB A key prerequisite to perform interferometry over long baselines is the 
preparation of a very narrow velocity distribution even beyond the ones of typical Bose- 
Einstein condensates which was already demonstrated for both species [32] 33] 1331 [35]. 
This can be reached by delta-kick cooling [3.6, 37]. The facility we want to employ for a 
test of the UFF is the VLB AT Teststand located at the new founded Hanover Institute 
for Technology (HITec) [38]. This device will provide two experimental chambers for 
the preparation of atomic ensembles with two independent source chambers for a max¬ 
imum flexibility in the choice of atomic species. A 10 m ultra-high vacuum-tube with 
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Table 1 . Comparison of choices for test masses A and B employed in existing 
and planned tests of the UFF parametrized for violation scenarios with respect to 
their effective charges Qab an d /g e +p+™) fp e +v—fp g +p~ n > f^e+p+n calculated 

according to m and m Nuclide data is used from m and for Ti a natural occurrence 
of isotopes is assumed [22]. 


A 

B 

Ref. 

aq2, b 

■10 4 

AQk 2 ,n 

•10 4 

A/_ n 

•10 2 

A f +n 
•10 4 

A/_ n 

•10 5 

A/ +n 

•10 4 

9 Be 

Ti 

0 

-15.46 

-71.20 

1.48 

-4.16 

-0.24 

-16.24 

Cu 

238 pj 

DP 

-19.09 

-28.62 

-7.08 

-8.31 

-89.89 

-2.38 

6 Li 

7 Li 

ra 

0.79 

-10.07 

-7.26 

7.79 

-72.05 

5.82 

85 Rb 

87 Rb 

0123125] 

0.84 

-0.79 

-1.01 

1.81 

1.04 

1.67 

S7 Sr 

88 Sr 

m 

0.42 

-0.39 

-0.49 

2.04 

10.81 

1.85 

39 K 

87 Rb 

HU 

-6.69 

-23.69 

-6.31 

1.90 

-62.30 

0.64 

87 Rb 

170 Yb 

[This work] 

-12.87 

-13.92 

-1.36 

-8.64 

86.00 

-5.46 


a magnetically shielded region of approximately 9 m forms the baseline for an extended 
free fall. Since operation of the equivalence principle test only occurs in the magneti¬ 
cally shielded region we anticipate a free fall time of 1 s and up to 2.6 s if the atoms are 
launched. Assuming a measurement with 1 • 10 5 ytterbium atoms and 2 ■ 10 5 rubidium 
atoms produced in 10 s, this leads to a shot noise limited performance of 1.6 TO -10 Hz~ 1//2 
and 6.5 • 10~ 12 Hz -1 / 2 in the Eotvos ratio respectively. The second value relies on higher 
order beam splitters, as explained in chapter [5] 


4. Concept for a dual species source of rubidium and ytterbium 


Mixtures of rubidium and ytterbium have been studied before in various experiments 
[39] m but were not yet used for precision interferometry. The construction of a dual 
species source capable of supporting an EEP test experiment faces a variety of challenges 
which are studied in the first phase of the experiment described in this work. A source 
has to fulfill the following characteristics: 


• The clouds have to be able to be cooled down to quantum degeneracy to fully 
exploit the long time of free fall achievable in the used infrastructure. Although 
this is relaxed by employing so called delta kick cooling, the efficiency of this process 
is strongly dependent on the initial temperature. 


The initial collocation has to be very well known and controlled. To a certain degree 


this excludes isotope combination which arc immiscible as discussed in chapter 4.4 


• The initial velocity distribution of the two species has to be matched to a high 
degree to allow for differential suppression of systematic effects, like wave front 
curvature or residual rotations. 
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(a) Mach Zehnder geometry 



(b) Setup 


Figure 1 . Mode of operation in Mach-Zehnder configuration and sketch of the 
experimental setup. Shown in |1 (b)| is an operation in drop configuration. 


• To achieve the target performance, 1 • 10 5 ytterbium atoms and 2 • 10 5 rubidium 
atoms have to be brought to degeneracy in less than 10 s. If this performance is 
not reached, it will increase the time needed for integration, but is not prohibitive 
to the overall experiment. 

4-1. MOT Operation 

Rubidium has two stable isotopes with mass numbers 87 and 85, both are bosonic and 
can be brought to degeneracy with common methods [32], 33]. Since both are also 
natural abundant and can be cooled similar well by standard laser cooling techniques, 
the specific decision for a rubidium species will be taken based on the miscibility with 
the ytterbium isotopes. The widely spread method for the preparation of rubidium 
ensembles is laser-cooling on the 5 2 Si/2-5 2 P 3 /2 transition with a subsequent optical 
molasses step for achieving sub-Doppler temperatures down to approximately 2 pK. 
With a combination of a multi-layer atom chip allowing for an efficient transfer of laser 
cooled atoms to a magnetic trap and a 2D + -MOT, quantum degenerated ensembles with 
4 • 10 5 rubidium atoms were produced in 1.6 s [41 ] . 

With in total five bosonic and two fermionic stable isotopes that have all been brought to 
quantum degeneracy before [541135] . ytterbium offers a variety of choices for test masses 


















Testing the UFF with Rh- Yb in VLBAI 


7 


as seen in table [2] The bosonic isotopes have no hyperfine splitting and therefore a very 
low magnetic sensitivity compared to rubidium for example [42]. While this is beneficial 
to counteract systematic effects, the missing possibility to drive Raman-transitions 
between the hyperfine states is limiting the implementation scenarios. Ytterbium, an 
alkaline earth like element, offers the possibility to perform narrow-line cooling on the 
inter-combination transition l So~ 3 Pi with a Doppler-temperature of 7R = 4.4/iK. Due 
to a low vapor pressure one has to face the challenge to pre-cool the hot source for 
efficient MOT operation. The common method is the use of a Zeeman-slower with a 
transversal cooling stage at the singlet transition 1 Sq- 1 Pi |43j . Another comparably new 
option is the use of 2D-MOT at same transition [44]. Experimentally loading rates of 
6-10' 1/4 Yb atoms per second have been achieved by both methods. The 2D-MOT 
seems preferable over the Zeeman-slower setup in terms of vacuum quality in the main 
chamber due to the use of differential pumping stages and offers higher scalability with 
available laser power at 398.9 nm. 

Table 2. Stable isotopes of ytterbium and their relative natural abundance [43] in %, 
character of spin-statistic, intra-species scattering length jj43j, inter-species scattering 
length with 87 Rb in ao m , as well as isotope-shift relative to 174 Yb of the relevant 
cooling transitions in MHz. 


Isotope 

Abund. 

Spin st. 

(lYb/Yb 

a Yb/Rb 

J 

1 s 0 - 3 P 1 

1 S 0 - 1 P 1 

168 yb 

0.13 

boson 

252 ±3 

39.2 ± 1.6 


3655 

1887.4 

170 Yb 

3.05 

boson 

64 ±2 

-11.5 ±2.5 


2287 

1192.4 

171 Yb 

14.3 

fermion 

-2.8 ±3.6 

58.9 ±4.4 

(1/2-1/2) 

-2132 

1153.7 






(1/2-3/2) 

3805 

832.4 

172 Yb 

21.9 

boson 

-599 ± 64 

-161 ± 11 


1000 

1887.4 

173 Yb 

16.1 

fermion 

199 ± 2 

626 ± 88 

(S/2-5/2) 

2312 

-253.4 






(S/2-7/2) 

-2386 

588 






(5/2-3/2J 

3806 

516 

174 Yb 

31.8 

boson 

105 ± 2 

880 ± 120 


0 

0 

176 Yb 

12.7 

boson 

-24 ±4 

216.8 ±4.7 


-955 

-509.3 


4-2. Trapping and evaporation 

Since we aim for a combined trap of both species, magnetic traps are not an option for the 
magnetically not trappable ytterbium. As a result a far detuned optical dipole-trap in 
the mid-infrared will be used as a common trap. Figure [2] shows the scalar polarisability 
at a certain wavelength with respect to the inter-combination MOT for ytterbium. The 
differential polarisability shows mainly two remarkable results: Ytterbium is not trapped 
at 1 /mi and there is a zero-crossing close to 1.5 /mi, that would potentially allow for AC- 
Stark shift compensated dipole-trap. A more conservative and less demanding solution 
would be the use of a dipole-trap beyond the zero-crossing for example at 1960 nm. To 
compensate AC-Stark shift dispersion over the cloud, which would be large due to the 
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(b) Scalar polarizability 3 Pi 
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(d) Differential AC-Stark shift 


Figure 2. Scalar polarisability and effective AC-Stark shift. The upper curves 2(a) 


and |2(b)j show the laser wavelength dependent scalar polarisability of the states in the 
transition used for the intercombination line cooling. The lower curves show in 2(c) the 
differential polarisability and in |2(d)| the resulting differential AC Stark shift imposed 
on the intercombination line by a 1960 nm ODT with 100 W, a 50 /im waist and using 
an additional 8.84 mW dressing beam with 1 GHz blue detuned to the transition. 


narrow linewidth of the transition a low-intensity blue detuned compensation beam can 
be used [48] with a detuning of A comp . =2ir ■ l GHz and a power of I comp . = 8.84 mW. 
The Bose-Einstein condensation in a single beam dipole-trap at this wavelength for 8 'Rb 
was already shown in a weak hybrid trap configuration in |49j. Therefore, a 1960 nm 
trap appears to be an ideal solution and lasers with output powers up to 100 W are 
available. 

4-3. Dual species loading sequence 

The cycle time of the experiment will be limited by smaller loading rates of the 
ytterbium, even with the use of a 2D+-MOT and the expected increase in flux, due 
to the use of higher laser power. In addition the 1 Sq- 1 Pi transition cannot be driven 
together with the rubidium cooling transition 5 2 Si/2-5 2 P^/2, since the ionization energy 
of the upper state of rubidium is 2.59 eV that corresponds to 478.7nm. Therefore the 
dual species sequence will first completely undergo the loading steps for cooling and 
trapping ytterbium into the dipole trap before we start the fast loading of the rubidium 
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MOT. To avoid losses due to collisions at this stage of the experiment it is possible to 
shift the center of the rubidium MOT against the dipole trap via adjusting the magnetic 
field gradient before both isotopes are co-located inside the dipole trap. 

4-4- Species miscibility and dynamical evolution 

This ability to cool non-magnetic ytterbium isotopes to quantum degeneracy inside the 
2 //m dipole trap via evaporation without additional effort is a key motivation for our 
choice. Fermionic isotopes are not considered in this study since degenerate Fermi gases 
are large and expand with higher rates than BECs, which is an important parameter 
for long baseline interferometry. They might nevertheless be interesting for future tests 
and the device is designed to keep this option open. As table 2 shows, we are left with 
five bosonic isotopes where two of them, 1,2 Yb and 176 Yb, have negative intra-species 
scattering length. They would require a more complex experimental design including the 
manipulation of an optical Feshbach resonance to reach degeneracy. 174 Yb is the most 
abundant isotope which was already condensed [35] . Nevertheless, due to the repulsive 
collisions to 8 'Rb (inter-species scattering lengths of (880 ± 120) %), a binary mixture 
will not be stable due to three-body losses. For all the reasons stated above, we focus 
our investigations on 168 Yb, 1,0 Yb and possible mixtures with 8 'Rb. Unfortunately, 
168 Yb and 170 Yb are the least abundant isotopes making loading rates significantly 
low which constrains the cycling rate in the order of tens of seconds unless they are 
enriched. The 168 Yb - 8 'Rb mixture features an inter-species positive scattering length 
of 39.2 ± 1.6 Oo meaning that this Yb isotope can be sympathetically cooled by 8 'Rb 
atoms. As shown in our systematics study in section [5j the separation between the 
two components of a binary mixture has a dramatic effect on the performance of the 
UFF test. Therefore, quantum miscibility cannot be neglected in this density regime. 
Indeed, if the interspecies repulsion exceeds the miscibility threshold [50j, the two atomic 
clouds spatially separate to minimize the interaction energy. This immiscible state is a 
hindrance for optimising the overlap of the centre of mass of the two wave packets fed 
into the interferometer for comparison. This makes it necessary to carefully check for the 
proposed isotopes if they can be prepared in overlapping pairs of spherical symmetry. 
We therefore solve a system of 3D-coupled Gross-Pitaevskii equations describing the 
ground state of the mixture {5Tj . The results of these simulations are shown in figure [3] 
The calculations confirm the miscibility of 8 ‘ Rb with the two Yb isotopes considered 
making it a suitable candidate for an UFF test. In contrast, the combination of 168 Yb 
with 1,0 Yb builds up a symmetric shell structure. These binary states numerically found 
are susceptible to and deformable by external fields (magnetic forces, gravitational sag, 
etc.) present in the science chamber. Therefore, this mixture is not considered for 
dynamics and systematics. 

In order to reduce systematic errors of the atom interferometric comparison and allow for 
an extended interrogation time, it is crucial to reduce the size of the atomic samples. In 
the proposed facility, few seconds of free fall or launch time are used to reach the target 
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(a) 168 Yb- 87 Rb (b) 170 Yb- 87 Rb 

Figure 3. Density plots of the ground states of the 1,0 Yb/ 168 Yb and s 'Rb mixtures. 
For each pair mixture, the wave functions are computed solving the Gross-Pitaevkii 
equation in 3D including the intra-species interactions of the two isotopes and the 
inter-species one with 8 'Rb. The magnitudes of these interactions are the same shown 
in table [2j We assume that each mixture is confined by the same external trap with 
frequencies solely differing due to the mass difference. The trapping frequencies are 
2n ■ 88 Hz for Rb and 2i r • 67 Hz for Yb. In both cases, a symmetric mixture ground 
state is found illustrating the miscibility of the two pairs without further tuning of 
external optical or magnetic parameters (Feshbach for example). 


accuracy of the UFF test. It is clear that thermal ensembles would reach very large sizes 
at these time scales. This motivates the use of degenerate matter waves characterized 
by a slow expansion. The state-of-the art in slowing down the expansion of BECs 
improved dramatically with the use of delta-kick cooling (DKC) techniques [36] [12]. In 
recent experiments with a comparable baseline [37], it was experimentally demonstrated 
that the expansion energy of a degenerate 8 'Rb ensemble could be restricted to only 
few tens of pK in 2D. We anticipated such records when proposing space missions with 
more than 10 s of free evolution time m of a mixture of a 8 'Rb / 85 Rb condensates. 
The DKC manipulation [52] consists in collimating matter waves by suddenly reducing 
the frequency of the initial trap holding the atoms and cutting it when all atoms reach 
the turning points of the trap walls (at t p /4, where t p is the trap period). The same result 
is expected by re-pulsing the initial trap after switching it off for some free expansion 
time. A substantial part of the atoms kinetic energy is absorbed by this process leading 
to a slowed expansion. The analogy with light beams collimation often led to label this 
manipulation as an atomic lens. We anticipate the use of a double lens to match the 
expansion rates of 8 'Rb an 1,0 Yb. This match is mandatory to mitigate errors related 
to residual wave front curvatures and relaxes the requirements on the initial collimation 
and retro reflection mirror planarity. 

4-5. Interferometer sequence 

As described earlier, performing an UFF-test is equivalent to a simultaneous 
measurement of the gravitational acceleration g^B acting on the two test masses. To 
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perform this measurement with atoms a sequence of light pulses has to be applied to 
interrogate them with respect to a common reference mirror which acts as a phase front 
reference. The most prominent configuration for inertial sensitive atom interferometry 
is the Mach-Zehnder-type 7 t/2 — i r — 7 t/ 2 sequence with a time T of free evolution in 
between each of the pulses. Two different modes of operation can be distinguished: (i) 
dropping atoms from a source on the top of the device and (ii) launching atoms onto 
a parabolic trajectory from a source at the bottom of the device. While the first mode 
is characterized by a good control over the initial conditions at free evolution times 
of 2 T = 1 — 1.3 s at a baseline of roughly 9 m, the second one offers the perspective 
to increase the overall length of the interferometer up to 2 T = 2.6 s. Launching over 
approximately 10 m was already demonstrated for rubidium in an accelerated optical 
lattice by coherently transferring a large number of photons at a decent efficiency |12j 
and appears also realizable for ytterbium with similar parameters. Nevertheless, this 
fountain mode requires a well controlled launching velocity of both test masses. 

4-6. Beam splitting and match of scaling factor. 

A major limitation for inertial measurements with atom interferometers is seismic noise 
which scales similar to the acceleration signal with T 2 and thus limits the maximum time 
of interferometry where the signal to noise ratio is still improving. When using a common 
mirror for a differential measurement, as planned for this experiment, the seismic noise 
for both interferometers is common and thus suppressed in the difference signal [53l 53]. 
To fully benefit from the non magnetic properties of the ytterbium 1 S'o state and allow 
for higher order beam splitting we plan to use Bragg type beam splitters, coupling 
momentum states of the respective ground states. The used off resonant transitions 
are the 1 Sq- 1 P 1 transition for ytterbium at 399 nm and the 5 2 S'i/ 2-5 2 P3/2 transition for 
rubidium at 780 nm. The suppression factor depends on the match of the scaling factor 
kT 2 , with the effective wave vectors k , and of the sensitivity function which is itself 
dependent on the timing of the interferometer pulse sequence. The basic approach 
is to match the scaling factors by tuning the interferometry time T for each species 
individually [52]. This will lead to a small difference in the frequency response of the 
two interferometers and will not properly suppress contributions scaling differently with 
T but allows for a simple data analysis scheme. 

In the case of mismatched effective wave vectors and same pulse timing, the phase 
frequency response is similar between the two species but rescaled according to the 
appropriate wave vector. As long as the resulting phase noise is smaller than 1 rad the 
phase information can still be fully recovered by weighting the results with the wave 
vector ratio. An analysis of this case can be found in [53]. Even in the case of noise 
above tt most of the information can be recovered at the cost of signal to noise ratio. In 
the case of higher common noise contributions the resulting 2n ambiguity can be fully 
resolved by operating an additional classical sensor [15]. Another option is to adapt 
the model used for data interpretation and recover at least some level of suppression by 
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5. Requirements and error budget 

This chapter summarizes the requirements on experimental and environmental 
parameters to restrict statistical and systematic errors. These requirements are partly 
relaxed compared to single species gravimetry measurements [56j [57], because the 
simultaneous operation of the dual atom interferometer and certain parameters choices 
allow to engineer suppression ratios for inertial phase shifts and inhomogeneities in the 
beam splitting wave fronts. A detailed derivation and discussion of error terms for an 
UFF-test with 8 'Rb / 85 Rb in the 10 m tower in Stanford was reported in [3SJ and the 
error budget for a satellite based test can be found in [16], 159]. This paper utilizes the 
same approaches for error assessment and thus focuses on the results. 

We consider three different scenarios. In the near future, atoms will be dropped from 
the top chamber, and the scaling factors k^T^ = k Y bTy b will be matched. In this case 
of matched scaling factors, correlation between the two atom interferometers will then 
allow to extract the differential phase corresponding to the differential acceleration via 
ellipse fitting [_53J, iRD]. The next intermediate step is to use the same free evolution time 

= Tyb which mitigates bias terms ~ kT 3 , ~ kT l but requires a more complex 
read out scheme. Since the scale factors differ now, the correlated signal will not 
form an ellipse. Restricting phase excursion to below 2tt still allows the extraction 
of the differential phase via fitting the Lissajous figure [53]. However, the expected 
vibration noise level is above 2i r. As mentioned earlier this ambiguity may be lifted via 
correlation with classical sensor mounted in close proximity to the retro reflection mirror 
as demonstrated for an atom interferometer on a plane [T^j or by adapting the phase 
extraction algorithms. Finally, the advanced scenario considers launched atoms from the 
bottom chamber and increased momentum transfers by the beam splitters. A lattice 
launching technique inside a 10 m fountain ra and high momentum transfer beam 
splitters [61, [62] which meet the requirements of this paper were already successfully 
implemented by other experiments. Requirements for systematics are summed up in 
table [3] and the resulting uncertainties in table [4] Statistical fluctuations in these 
parameters are allowed up to the levels reported in table [5] which implies the errors 
in table IB 

To engineer a high common mode rejection ratio, the center of mass positions, 
center of mass velocities, size and expansion ratios of the two atomic species have to be 
matched. Coupled to gravity gradients and rotations, position and velocity differences in 
the center of mass positions cause spurious phase shifts in the differential signal. Using 
trapping frequencies of 27T • 500 Hz implies a gravitational sag of 1 pm which will need to 
be characterized to 1 % in the advanced scenario. Due to the lattice launch, we expect 
a differential velocity of 31 pm/s. The corresponding biases will be subtracted from the 
signal which imposes the requirement of knowing the gravity gradient to 0.1%. This 
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Table 3. Requirements to reach the stated uncertainty in g in different configurations. 
1) Assuming T zz = -2.91 • 10 -6 s~ 2 , \T XX \ = \T yy \ = \T zz \/2, and T zzz = 5 Z T ZZ . 2) 
Assuming a counter rotation 12 c with quality (£l c —Q y )/n y , Q y = 57.5 grad/s. 3) Waist 
after collimation lens wq, ideal distance lens - fiber 40 cm, assumed distance 40cm-df. 


Error source 

Initial 

Intermediate 

Advanced 

Free evolution time 

500 ms, 

500 ms, 500 ms 

1300 ms, 


505.7 ms 


1300 ms 

Effective wave vectors k Rb , 

2-47t/( 780 nni 

l), 47r/(399nm) 

8-47r/(780nm), 

k Yb 



4-47t/( 399 nm) 

- Relative uncertainties 

±10 8 in individual k, 


±10 14 in the ratio ( ky-b ~ km 

)/{kyb ± k R b) 

Rabi frequency match 


1% 


Common velocity v 0 

-3 m/s (±30 pm/s) 

12.75 m/s 




(±0.1275 prn/s) 

Differential position 

1 pm 

1 pm±0.1 pm 

1 pm±10 nm 

Differential velocity 

1 pm/s 

0.1 pm/s 

31 pm±10 nm/s 

Grav. gradients - T zz ^ 

3 • HR 6 (±1.5 • 

3 • 10" 6 (±3 • 

3 • 10~ 6 (±1 • 


10” 10 ) s -2 

10- 7 )s' 2 

10 -9 )s” 2 

Grav. grad. - T zzz 

< 5-10 -9 (±3.3- 

< 5-10 -9 (±3.3- 

< 5 • 10“ 9 


10 -11 ) nr 1 s~ 2 

10 -11 ) m _1 s~ 2 

(±10~ n ) m _1 s~ 2 

Acceleration compensa- 

10 -7 ms -2 

5 • 10 -8 ms -2 

10 -8 ms -2 

tion via frequency scan 




to 




Counter rotation quality 
2) 


2% 


Rotation 


1 prad/s 


Rotation Yt z 


44.4 prad/s 


Inside shield 




- Magnetic offset field B 0 

25 mG 

25 mG 

25 mG 

- Magnetic field gradient 

50 pG/rn 

15 pG/m 

1.5 pG/rn 

SB 




Inside source chamber 




after release Bo, AB 


< 1G, < 0.1 G/m 

Collimation df, Wo 

100pm, 2.3 cm 

100 pm, 3 cm 

Mirror quality 


A/20 


Initial sample radius 


300 pm 


Effective temperature 

25.6 nK(±l %), 

2.5 nK(±l %), 

256pK(±l %), 


50 nK(±l %) 

5nK(±l%) 

500pK(±l %) 

Atom numbers 

2-10 5 (±1%), 10 5 (±1%) 

Scattering lengths 

(100.4±0.1)a 0 [55], (64±2) a 0 , (11.5±2.5)a 0 

Beam splitting accuracy 

0.01 

0.001 

0.001 
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Table 4. Contributions of the different error sources to the uncertainty in g in different 
configurations. 1) Requires back correction via knowledge of g, T zz , T zzz , and Q y 


Error source 

Initial 

Intermediate 

Advanced 

u, ? 

in 10- 12 

o 

1 

CO 

10~ 14 

Gravity gradient + posi- 

0.3 

0.3 

0.3 

tion overlap 




Gravity gradient + veloc- 

0.15 

0.15 

0.4 

ity overlap 




Gravity gradient + g, v 0 

0.15 

0.15 

0.15 

Coriolis x 

0.23 

0.23 

0.23 

Coriolis y 

0.2 

0.2 

0.2 

Other terms d 

1 

1 

1 

Magnetic fields 

0.3 

1 

1 

Wave fronts 

5.1 

5.2 

5.7 

Mean field 

1.3 

3.6 

3.9 

Sum 

5.7 

6.7 

7.4 


Table 5. Requirements on noise sources for the dual species atom interferometers 
in different configurations. All contributions are expected to be uncorrelated. The 
requirements were set to reach the shot noise limit. Where appropriate values are 
given as a requirement for a single measurement cycle. (1) Assuming correlation 
with an additional classical seismometer or advanced data fitting eliminating the 2ir 
ambiguity. 


Noise source 

Near / intermediate 


Advanced 

Shot noise 

See tab. [5 

for N, k, and T. 

Beam splitter 

1 kHz Lor< 

mtzian linewidth 

Linear vibrations 

10 -6 ms -2 Hz _1/2 


10 -6 ms -2 Hz _1/2 d) 

Starting velocity 

a v < 0.3nmi/s 


a v < 3.8 prn/s 

Overlap 

<7Ar < 10 pm, 


OAr <0.3 pm, 


a Av < 10 pm/s 


(j Azj < 0.3 pm/s 

Magnetic fields 

(Jsb < 0.5mG/m 


a$B < 45 pG/m 

Wave fronts 

crdf — °Az = 100 pm, jitter telescope & mirror position 1mm 


in z-direction (g) 

Mean field 

5 % jitter in beam splitting 

1 % jitter in beam splitting 


ratio, 20 % in atom numbers 

ratio, 20 % in atom numbers 

Cycle times 

11s 


12.6s 


will be measured with the apparatus itself in a gradiometer operation mode. Existing 
gradiometer experiments reached noise floor of down to 3 • 10~ 8 s~ 2 Hz -1 / 2 [63j f28] . 
Furthermore, a counter rotation of the retro reflection mirror will reduce the bias due 
to the earth’s rotation na. Additional errors occur if the atoms map different parts of 
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Table 6. Resulting noise contributions following tab. [5] All contributions are expected 
to be uncorrelated. The requirements were set to reach the shot noise limit. All values 
are given as the noise of a single measurement. 


Noise source 

Near / intermediate 

Advanced 


in 10“ 10 m/s 2 

in 10“ n m/s 2 

Shot noise 

4.8 

1.8 

Beam splitter 

2.8 

1 

Linear vibrations 

2.8 

1.8 

Overlap 

1 

0.3 

Starting velocity 

0.1 

0.03 

Magnetic helds 

0.3 

0.3 

Wave fronts 

0.12 / < 0.01 

< 0.01 

Mean held 

0.6 

0.4 

Sum 

6.3 

2.8 

- after 24 h 

- 7.1T0 -2 

- 3.4-10" 2 


the beam splitter wave fronts to which imperfect collimation or the finite quality of the 
retro reflection mirror cause inhomogeneities. Commercially available mirrors are rated 
up to A/20 (peak to valley) [64] which puts requirements onto the maximum allowable 
expansion rates. Demonstrated perfomances of lensing 8 'Rb atoms to InK in 3D [ 36] . 
and to 50 pK in 2D [37] are sufficient for the experiment. 

Additional sources for errors are magnetic holds inducing a second order Zeeman shift in 
the 8 'Rb interferometer and the scattering properties of the individual ensembles and the 
mixture. Suppression of magnetic stray holds with residual rms deviations of ~0.8mG 
inside a three layer 8.8m /i-metal shield were demonstrated [S3]. Therefore, additional 
calibration might be necessary to characterize the magnetic helds to the required level. 

6. Conclusion and outlook 

We presented a novel experimental scheme to test the EEP with two different atomic 
species, namely ytterbium and rubidium which is in the progress of being set up in 
Hanover in the new infrastructure of the Hanover Institute of Technology. Using this 
particular test pair for precision inertial sensing with atom interferometry imposes some 
challenges which are discussed in this letter together with appropriate specihc solutions. 
Based on the knowledge of this kind of measurement we provide an assessment of the 
expected performances of the experiment and of the major systematic effects. They 
should allow to test the Eotvos parameter at a level of 7 • 10 -13 in the next few 
years. The work described in this letter is the first step in a complete investigation 
of inertial sensing with an alkaline earth like element as ytterbium. In the framework of 
the collaborative research center geo-Q we will investigate possible applications of this 
technology for geodesy and further ways to improve ground based EEP tests beyond the 
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level of tests with devices employing classical test masses. We expect this work to have 
a major influence on to the field of fundamental sciences by giving new limits to possible 
violation scenarios. Moreover, the possibility to investigate interferometric techniques on 
long time scales with a high repetition rate will benefit atom interferometry experiments 
in micro-gravity environment or space platforms. 
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